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Abstract 

Three versions of the quark-meson coupling (QMC) model are applied to 
describe properties of nuclear matter and finite nuclei. The models differ 
in the treatment of the bag constant and in terms of nonlinear scalar self- 
interactions. In two versions of the model the bag constant is held fixed at its 
free space value whereas in the third model it depends on the density of the 
nuclear environment. As a consequence opposite predictions for the medium 
modifications of the internal nucleon structure arise. After calibrating the 
model parameters at equilibrium nuclear matter density, binding energies, 
charge radii, single-particle spectra and density distributions of spherical nu- 
clei are analyzed and compared with QHD calculations. For the models which 
predict a decreasing size of the nucleon in the nuclear environment, unrealistic 
features of the nuclear shapes arise. 



PACS number(s): 12.39.Ba, 21.60.-n, 21.10.-k, 24.85. +p 
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I. INTRODUCTION 



Although quantum chromo dynamics is believed to be the fundamental theory of strong 
interactions, low- and medium-energy nuclear phenomenology is successfully described in 
terms of hadronic degrees of freedom. Theoretical challenges arise in phenomena which reveal 
the quark structure of hadrons. Building models which connect observed nuclear phenomena 
and the underlying physics of strong interactions has become one of the fundamental goals 
in nuclear physics. Such models are necessarily crude and contain a variety of unknown 
parameters since the study of the nuclear many-body problem on the fundamental level is 
intractable. However, it is important that the new models respect established results which 
are successfully described in the hadronic framework. 

The quark- meson coupling (QMC) model proposed by Guichon JTJ, provides a simple 
framework to incorporate quark degrees of freedom in the study of nuclear many-body 
systems. In the QMC model nucleons arise as non-overlapping MIT bags interacting through 
meson mean fields. The model has been applied to a variety of problems in nuclear physics. 
It was shown that it describes the saturation properties of nuclear matter [0-0], and that it 
gives a fair description of the bulk properties of finite nuclei p|-I5|. The model was extended 



to include hyperons [L4|] and applied to studies of hyper-nuclei [p~5j,|T€ 



Although it provides a simple and intuitive framework to describe the basic features of 
nuclear systems in terms of quark degrees of freedom, the QMC model has a serious short- 
coming. It predicts much smaller scalar and vector potentials than obtained in successful 
hadronic models @J^,§]. As a consequence the nucleon mass is too high and the spin-orbit 
force is too weak to explain spin-orbit splittings in finite nuclei PHTT . 



The QMC model can be significantly improved by introducing the concept of a density 
dependent bag constant |5||| . The modified quark- meson coupling (MQMC) produces large 
scalar and vector potentials under the condition that the value of the bag constant in the 
nuclear environment significantly drops below its free-space value. 

In the present work we study three different versions of the QMC model. The first model 
is the original QMC model which we denote by QMCI. The second model denoted by QMCII 
is an extension which includes cubic and quartic self interactions for the scalar meson. As 
our third model we adopt model MQMC A from Ref. [[?[], in which the bag constant is a 
function of the scalar field. 

The original idea of the QMC model was to calibrate the model in free space such that 
the nucleon mass is reproduced and then to extrapolate to many-nucleon systems. However, 
it is not possible to account for all the necessary degrees of freedom. Most importantly, the 
scalar field which describes the mid-range part of the nucleon-nucleon interaction cannot be 
described in the framework of the simple bag model, but it must be included. The additional 



cubic and quartic scalar self-interactions of the QMCII model were introduced in Refs. [11 



to model the density dependence of the scalar mass. The corresponding couplings cannot 



be determined on a fundamental level pi| . In our approach we choose these parameters to 
reproduce the properties of nuclear matter near equilibrium that are known to be character- 
istic of the observed bulk and single-particle properties of nuclei. This procedure guarantees 
the reproduction of large scalar and vector potentials. 

In the MQMC model new parameters arise from the unknown density dependence of the 
bag constant. A refined MQMC model can be accurately calibrated to produce the empirical 
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saturation properties of nuclear matter M and it provides a good description of the bulk 



properties of finite nuclei [|13 . Most importantly, the model can be calibrated to reproduce 



small values for the effective nucleon mass which leads to a realistic description of spin-orbit 



splittings |13 



One of the basic motivations of applying quark models to nuclear systems is the hope to 
describe medium modifications of the internal quark-substructure of the nucleon. An impor- 
tant phenomenological quantity is the nucleon size which is represented in the framework of 



the QMC model by the bag radius p7| . Due to the increase of the bag constant the MQMC 
model predicts significantly swollen nucleons in the nuclear environment In contrast, 

for the QMCI and II model the bag radius decreases with increasing density P,p|,|9|,[TTH 

Our main goal is to study properties of nuclear matter and finite nuclei. We investigate 
whether the different versions of the QMC model lead to results which are consistent with 
established hadronic phenomenology. At present such an analysis might provide the only 
testing ground for the reliability of the strikingly different predictions for the changes of the 
internal nucleon structure. 

A well established framework for relativistic hadronic models is provided by quantum 
hadrodynamics (QHD) ||18|| . Numerous calculations have established that relativistic mean- 
field models based on QHD lead to a realistic description of nuclear matter and of the bulk 
properties of finite nuclei throughout the Periodic Table [18-25]. A very compelling feature 
of QHD has been the reproduction of spin-orbit splittings in finite nuclei. For comparison we 
employ a QHD model that includes quartic and cubic scalar self-interactions. We calibrate 
the model parameters so that QMCII, MQMC and QHD lead to the same nuclear matter 
properties at equilibrium. The QMCI does not provide a sufficient number of parameters to 
fit all the desired nuclear matter properties. Particularly, a rather high value for the effective 
nucleon mass at equilibrium is a prediction in this model 0. 

In our previous work 0,[y| we demonstrated that the QMC models are equivalent to a 
QHD type mean field model with a general nonlinear scalar potential and a coupling to the 
gradients of the scalar field. One of the key observations in the success of hadronic models is 
that nonlinear scalar self-interactions, i.e. a nonlinear potential, must be included |l9|-27 



This potential is constrained by nuclear observables and we check if the QMC models predict 
the typical size and form of the potential. 

In nuclear matter we find that the models which can be accurately calibrated, namely 
QMCII, MQMC and QHD, lead to essentially identical results. The MQMC and QMCII 
models reproduce large scalar and vector potentials and the typical size of the nonlinear 
potential. In contrast, the QMCI model leads to rather small scalar and vector mean fields. 
At the equilibrium point the nonlinear potential has only half the size of the value predicted 
by the other models. 

The results for the binding energies and spin-orbit splittings are similar. The QMCII, 
MQMC and QHD model, which reproduce the empirical properties of nuclear matter, lead 
to a realistic description of the experimental numbers; the only exception are the light nuclei. 
The surface energy in the QMCII model is too small and as a consequence the light systems 
are systematically overbound. The QMCI model underestimates the binding energies and 



gives only a poor reproduction of single-particle spectra and spin-orbit splittings PHTTf . The 
key observations is that the MQMC and QMCII model can reproduce sufficient small values 
for the effective mass, which is strongly correlated with the spin-orbit force in relativistic 
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mean-field models. 

The main difference between the models arises from the nonlinear coupling to the gradient 
of the scalar field. The coupling has no effect in nuclear matter and its size is a prediction 
of the underlying bag model. It mainly effects the nuclear shapes leading to more diffuse 
surfaces in the QMC models. At normal nuclear matter density the gradient coupling in the 
QMCII model is 6-7 times bigger than in the MQMC model. From the point of view of an 
effective field theory a gradient coupling arises naturally as a subset of possible nonlinear 
meson-meson interactions f25fl . However, rather than attempting to compete with these 
more sophisticated versions of QHD, it is our goal to analyze if an approach which is based 
on a simple quark model can reproduce well established results of nuclear phenomenology. 
For comparison we employ a more conventional version of QHD which does not include the 
gradient coupling. We compensate the effect on the nuclear shapes by adjusting the scalar 
mass to reproduce the experimental value of the charge radius in 40 Ca. Although after the 
adjustment all the models predict nearly identical rms charge radii we find sizable differences 
in the predicted density profiles. The QMCI and QMCII model lead to very compact nuclei 
with small central densities and steep surface areas. This effect is more pronounced for the 
light nuclei. 

The outline of this paper is as follows: In Sec. II, we review the formalism for finite nuclei 
and nuclear matter. Section III contains a short summary of the QMC model and the rela- 
tions which determine the properties of the nucleon. We also briefly discuss the calibration 
procedure. In Sec. IV, we analyze nuclear matter. We concentrate on the predictions for 
the nonlinear potential. In Sec. V, we apply the models to finite nuclei. Sec. VI contains a 
short summary and our conclusions. 

II. THE QUARK-MESON COUPLING MODEL FOR NUCLEAR SYSTEMS 

To study the properties of finite nuclei and nuclear matter we use a relativistic mean- 
field model containing nucleons, neutral scalar (0) and vector fields (V^) and the isovector 
p meson field (b M ). For a realistic description of finite nuclei, the electromagnetic field (A M ) 
must also be included. We assume the nucleons obey the Dirac equation 

(it-?- ~T- P - i(l + r 3 ) / - M*) i, s {x) = . (1) 

The potentials (V P ,B U ,A U ) and the effective mass M* are functionals of the meson mean 
fields, their form depends on the underlying quark model. 

In the QMC model the quarks are described by the Dirac equation 

(i?-gr- \a q p T- jb - (1 + 3r 3 )~ 4 - K - 9 q J\) M*) = > ( 2 ) 

where m q is the current quark mass. The quark wave function is subject to the bag model 
boundary conditions at the surface of the bag. Because quarks and nucleons interact with 
the meson mean fields, Eqs. ([!]) and (|J) define a self-consistent scheme for the description of 
the nuclear system. In infinite nuclear matter (A u = 0) the meson mean fields are constant 
and the potentials are given by 
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V = ^V = g p b v . (3) 

The effective mass is an ordinary function of the scalar field, i. e. 

M* = M*(4>) . (4) 

For a finite system the solution of Eq. (|I|) and Eq. (0) is rather complicated due to the 
variation of the meson mean fields over the bag volume. In consequence, the quark wave 
function and the ground state of a bound nucleon are no longer spherically symmetric ||TU| . 



To make a numerical solution feasible it is necessary to calculate the quark properties by 
using some suitably averaged form for the meson mean fields. Here we follow the prescription 
of 0,^,|11|] and replace the meson mean fields on the quark level by their value at the center 



of the nucleon bag, i.e. we neglect the spatial variation of the mean fields over the bag 
volume. In this local density approximation the potentials in Eq. (|XJ) are simply obtained by 
the corresponding nuclear matter relations given by Eqs. @ and Eq. @. The corresponding 
relation for the electromagnetic field is 

A v = eA v . (5) 

If we restrict considerations to spherically symmetric nuclei only the Vq component of the 
neutral vector field and the neutral p meson field (denoted by bo) contribute. The ground 
state energy of a nucleus can be written as 



e n = J2 E i+J dvfyvtpr+um-hvvo) 

i=occ. V 



r + <v>\ (6) 



-i[(V6 ) 2 + ^]-^(VA ) 2 



where Ei are the eigenvalues of the Dirac equation Eq. (|IJ). 

To analyze a general class of models, we introduce a nonlinear scalar potential of the 
form 

U B {4>) = \ml<? + i^ 3 + iA0 4 . (7) 

The actual mean field configuration is obtained by extremization of the energy. This leads 
to the set of self-consistency equations 

(A-m 2 v )V = -g vP , (9) 

(A-m 2 p )b = -g p - P3 , (10) 
AA = -ep p . (11) 

The densities on the right-hand side are the nuclear densities calculated with the wave 
functions in Eq. ([!]): 
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Ps = E , (12) 

i=occ. 

p= E ftYVi, (is) 

i=occ. 

p 3 = E ^>3 T >; , (14) 

i=occ. 

P P =^ E ^(1+73)7°^ ■ (15) 

i=occ. 

In the limit of infinite symmetric nuclear matter Eqs. (|6|) simplifies to 

% = - 2 t dk k\k 2 + M* 2 ) 1/2 + 9vV p - \m\V 2 + UM . (16) 
For the time-like component of the vector field Eq. @ reduces to 

9,V = 4p , (17) 

whereas the scalar field is determined by the equivalent of Eq. (g): 

dU* dM* 2M* /- fe F k 2 

o ^ (feWF ' (18) 



n 2 



In symmetric matter the Fermi momentum of the nucleons is related to the conserved baryon 
density by 

'=5?*- (19) 

The details of the underlying quark substructure are entirely contained in the expression 
for the effective mass M*((f>). In the next section we will discuss the functional form of the 
effective mass in the framework of the QMC model. 



III. THE QUARK-MESON COUPLING MODEL 

In this section, we briefly summarize the relations which determine the nuclear equation 
of state in the quark-meson coupling model. For further details we refer the reader to 
Refs. 

In the QMC model the nucleon in the nuclear medium is described as a static, spherical 
MIT bag in which quarks couple to meson mean fields. 

The energy of a bag consisting of three quarks in the ground state can be expressed as 

E ba9 = 3^ - | + \*I?B . (20) 

where the parameter Z accounts for the zero point motion and B is the bag constant. The 
coupling of the quarks to the scalar field is inherent in the quantities Q q and x which are 
given by 
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* 



- Rm*\ 1 ^ 2 
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and where m* = — denotes the effective quark mass and m" is the current quark 
mass. In the following we choose = 5 MeV. 

To remove the spurious center-of-mass motion in the bag we follow Ref. || and adjust 
the parameter Z in Eq. fl2"0| ) to reproduce the experimental value of the nucleon mass in the 
vacuum, i.e. we take 

M*(<j)) = E bag . (22) 

For a fixed meson mean-field configuration the bag radius R is determined by the equi- 
librium condition for the nucleon bag in the medium 



In free space M can be fixed at its experimental value 939 MeV and the condition Eq. ([23]) 
to determine the parameters B = Bq and Z = Zq. For our choice, Rq = 0.8 fm, the result 
for Bq^ and Zq are 169.97 MeV and 3.295, respectively. 

In the original version of the QMC model the bag parameters B and Z were held 
fixed at their free space values B = Bo, Z = Zq. Formally, the bag constant B is associated 
with the QCD trace anomaly. In the nuclear environment it is expected to decrease with 
increasing density as argued in Ref. p8| . 

To account for this physics in the QMC approach Jin and Jennings |||6| proposed two 
models for the medium modification of the bag constant: a direct coupling model in which 
the bag constant is a function of the scalar field and a scaling model which relates the bag 
constant directly to the effective nucleon mass. The density dependence is then generated 
self-consistently in terms of these in-medium quantities. In a previous work [^] this approach 
was generalized and it was demonstrated that the resulting improved MQMC model can be 
accurately calibrated to reproduce the empirical properties of nuclear matter |7| and finite 
nuclei jOJ. 

For our purpose here we adopt the model MQMC A of Ref. J7| in which the bag constant 
depends on the scalar field only 



v 



Tq = { 1 ~ g *M m ) With F{0) = 1 ■ (24) 

We model the functional form of F by using a simple polynomial parametrization 

F(<f>) = 1 + a<P + f34> 2 . (25) 

The functional form Eq. ( p5| ) provides sufficient flexibility and it is not necessary to include 
further nonlinearities in the scalar potential. We therefore analyze the MQMC model with 

U s {4>) = . (26) 
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In addition we consider the original QMC model in two versions. In the first model, which 
we denote QMCI, we disregard the nonlinear terms in the scalar potential and use the same 
form as in Eq. (|26|). In the second model, which we denote QMCII, we include the cubic 
and quartic terms as given by Eq. (|7|). We analyze if these additional nonlinearities lead to 
an improvement of the model. 

We close this section with a brief description of the calibration procedure. The param- 
eters Bq and Z are fixed to reproduce the nucleon mass in the vacuum. In nuclear matter 
the MQMC model contains seven free parameters. The parametrization of the bag constant 
contains the parameter r\ and the three couplings (g B , a, (3); in addition values for the ratios 
gl/m s ,g v /m v ,g p /m p are needed. 

To determine the parameters we follow @ . For given values of rj and gj / m s we determine 
the values of (g B ,a, /3, g v /m v , g p /m p ) to reproduce the equilibrium properties of nuclear 
matter, which are taken to be: the equilibrium density and binding energy (p°, — e ), the 
nucleon effective mass at equilibrium (M*) the symmetry energy (a 4 ) and the compression 
modulus (K ). The set of equilibrium properties used here are listed in the first row of 
Table | (denoted by input). For more details concerning the calibration procedure we refer 
the reader to Ref. 0. 

The model QMCII contains the five parameters g^/m s , g v /m v , g p /m p , k and A which are 
determined to reproduce the same equilibrium properties as for the MQMC model. In 
the model QMCI only the parameters g%/m s , g v /m v , g p /m p are available. We chose to fix 
these three parameters such that the model reproduces the same density, binding energy 
and symmetry energy as in MQMC and QMCII. Correspondingly, the effective mass and 
compression modulus at equilibrium are a prediction. The values are quoted in the second 
row of Table |. We find M*JM = 0.80 and K = 280 MeV. 

For finite nuclei calculations values for the meson masses are needed. The mass of the 
scalar meson m s is determined to reproduce the charge radius in Ca as we will discuss in 
more detail in section V. The masses of the remaining mesons are fixed at their experimental 
values m v = 783 MeV and m p = 770 MeV. 



IV. PROPERTIES OF NUCLEAR MATTER 

Quark-meson coupling models are designed to describe both bulk properties of nuclear 
systems and medium modifications of the internal structure of the nucleon. It is impor- 
tant that the models reproduce established results of nuclear phenomenology before reliable 
predictions for changes of the quark substructure can be made. 

An important quantity is the nucleon size which is represented in the framework of the 



QMC model by the bag radius [|17]]. The density dependence of the bag radius can be seen 
in Fig. |l[ The opposite behavior of the prediction of the models is striking. The MQMC 
model leads to significantly swollen nucleons 0-0. At equilibrium the bag radius increases 
to roughly 25% of its free space value. At low and moderate densities we observe a very 
small dependence on the model parameters. In contrast, for the QMC models the bag radius 
decreases slightly with increasing density P,^|flpT| . The effect is more pronounced for the 
model QMCII. 

The medium dependence of the effective nucleon mass is of central importance in rela- 
tivistic nuclear phenomenology. The effective mass is shown in Fig. as a function of the 
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density. Also indicated is the QHD result based on a model which includes a nonlinear scalar 
potential of the form given by Eq. ([/]). The QHD parameters are determined by reproducing 
the equilibrium properties listed in Table |. The masses for MQMC, QMCII and QHD are 
nearly identical at low and moderate densities. In contrast, QMCI predicts a very high and 
slowly decreasing effective mass ||. 

The bag models cannot be extrapolated to arbitrary high densities. The solutions cease 
to exist when the point x = with Rm* = —3/2 in Eq. fl2"T| ) is reached. This corresponds 
to a maximal density p max ~ 4.92p° and p max « 1.23p° for the QMCI and QMCII models 
respectively. The rather small value for the QMCII model severely limits its applicability. 
The MQMC model leads to p max ~ 3.38p°. Here applications are limited by the large bag 
radii rather than by the maximal density. The individual bags start overlapping before p max 
is reached [0]. 

As discussed in Ref. there is a direct relation between the QMC models and QHD- 



type mean field models. The main difference between QMC and QHD is the functional form 
of the effective mass. In QMC it is a complicated function of the scalar field 

^QMC = KmM . (27) 

whereas in QHD it is linearly related to the scalar field 

M* HD = M — g $ . (28) 
This suggests a redefinition of the scalar field in QMC 0: 

gM4>) = M - M* MC (0) = M-(3^ - | + ~7tR 3 B) . (29) 
The coupling g is chosen to normalize the new field according to 

$(0) => + O(0 2 ) , 



and is given by 

dM*((f)) 



9o 



(30) 

<t>=o 



The scalar potential in Eq. ([?]) can now be expressed in terms of the new field 

u s m = u s (m) = \rn 2 s <t?m + ^ 3 ($) + . m 

The field redefinition Eq. (p9|) recasts the QMC model into a QHD type mean field model. 
Most importantly, it predicts a specific form for the nonlinear scalar potential. Fig. ||| 
indicates the predicted potentials as a function of the transformed scalar field given by 
Eq. (||). Below the saturation point (g $/M = 0.37) the curves for MQMC and QMCII are 
almost identical to the QHD result. For a given value of the scalar field the QMCI potential 
is considerably bigger. At the saturation point (g Q/M = 0.2), however, it is only half the 
size of the other models. 
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The agreement in Fig. ^ is somewhat misleading. Although there is good agreement 
between QHD, MQMC and QMCII the functional form of the potentials is different. This 
can be studied by expanding the potential in Eq. (RTJ) in terms of the scalar field 



U B {</>{$)) = \e 2 (g^) 2 + ^e 3 (go<l>) 3 + ^ 4 (<7o$) 4 + • • • • (32) 

Fig. [| indicates U s and the series in Eq. ( |32| ) truncated at second, third and fourth order. 
Part (a) shows the result for QMCI and QMCII. The potential is well represented by the 
fourth order polynomial in Eq. fl32|). The major contribution comes from the quadratic 
term. The very small deviation of the truncated fourth order series from the exact potential 
indicates that higher order corrections are negligible. At the saturation point we find that 
the neglected higher order terms give a correction of only 2%. The result for the MQMC 
model is indicated in part (b) of Fig. f|. Here the third and fourth order contributions are 
more important than for the QMC models. At the equilibrium point the neglected higher 
order contribution give a correction of roughly 3%. 

It is well known that nonlinear self-interactions of the form Eq. (|7|), or Eq. (0), must be 
included for a successful low-energy nuclear phenomenology [|19H 27|1. Our analysis demon- 
strates that once the QMCII and MQMC models are calibrated by using characteristic 
properties of nuclear matter they predict the typical size of these nonlinear self-interactions. 
Furthermore, all the models appear to be natural, meaning that the nonlinear potential is 
well represented by a low order polynomial^]. 

The details of the underlying quark structure are entirely contained in the coefficients of 
the series in Eq. (|3T|). Based on a nuclear matter analysis alone, where the MQMC and the 
QMCII model produce equivalent results, it is not possible to decide which prediction for 
the changes of the internal nucleon structure is more reliable. 



V. CONSEQUENCES FOR FINITE NUCLEI 

The field redefinition in Eq. (|29|) can also be applied in a finite system. Expressed in 
terms of the new field $ the contribution of the scalar field to the energy in Eq. (|6]) is given 

by 

E S = J cn/[i(V0) 2 + £/ s (0)] (33) 

= / dv[^) 2 H^) 2 + u s m , 

In addition to the nonlinear potential Eq. Q3"T| ) the transformation also induces a coupling 
to the gradients of the scalar field 

H{*) = 1 +h(*) = || = -^. (34) 



1 For a more complete discussion of naturalness in the QMC model, see Ref. 2S] 
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The coupling has no effect in nuclear matter calculations and is a prediction of the 
underlying bag model. From a modern point of view, our model contains a subset of possible 
nonlinear meson- meson couplings. In more sophisticated versions of QHD p3[ , inspired by 
concepts and methods of effective field theory, these terms and many others are considered. 
For comparison we employ a conventional version of QHD which contains the standard form 
of the nonlinear potential given by Eq. (0) and with = 0. 

The coupling is indicated in Fig. [5]. Relevant for applications to finite nuclei is the 
region below g $/M = 0.37 for QMCII and MQMC and below g $/M = 0.2 for QMCI 
which corresponds to the saturation point of nuclear matter. The MQMC model predicts 
very small values for the gradient coupling. The function h(Q) is identical for the two 
versions of the QMC model. Here the coupling decreases significantly. Near nuclear matter 
equilibrium the coupling is 6 to 7 times bigger than for the MQMC model. 

The gradient coupling has an impact on the bulk and single particle properties of a 
nucleus. For positive values it effectively decreases the coupling strength of the scalar density 
to the scalar field $ leading to significantly smaller mean fields. The most prominent effect 
are changes of the nuclear surface. This can be studied in Fig. ^| which indicates the charge 
density for 40 Ca. Charge densities and charge radii are calculated by convoluting the point 
proton density Eq. ( |To] ) with an empirical proton charge form factor [|30|. In part (a) the 
same mass for the scalar meson (m s = 500.8MeV) was used for all the models. Note that 
all the models, except QMCI, predict the same equilibrium properties of nuclear matter. 
The gradient coupling changes the surface drastically. The density becomes smaller in the 
interior region and the surface area is more diffuse. As can be expected from Fig. |5| the 
effect is more pronounced for QMCII. To compensate for this effect we adjust the scalar 
mass such that the models reproduce the experimental charge radius in 40 Ca. The resulting 
charge densities are indicated in part (b) of Fig. Also included is the experimental charge 
density |[31|| . The adjusted values for the scalar mass are an indication of the size of the 
gradient coupling. The masses had to be increased by 10 and 30% for the MQMC and 
QMCII models respectively. After the adjustment the curves for QHD and MQMC are 
almost identical. Due to the very large value of the scalar mass the QMCII model leads to 
very compact nuclear shapes. The density is still too small in the interior region and the 
surface area is now very steep. Also included in Fig. || is the result for QMCI. Here the mean 
fields are much smaller than in the other models and the impact of the gradient coupling is 
weaker than for QMCII. In order to reproduce the charge radius of 40 Ca the scalar mass has 
to be decreased. 

The effect of the gradient coupling on the nuclear shapes depends strongly on the mass 
number of the system. Discrepancies in the predictions of the models are much smaller for 
the heavier nuclei. The quality of the reproduced nuclear shapes for the lighter nuclei can be 
studied in Figs. Fig. [7| indicates the isotope shift in calcium. Shown is the difference of 
the 48 Ca and 40 Ca charge density. The QMCI and QMCII model cannot provide a realistic 
description of the experimental curve ||32|| . In Fig. | we show the charge density of 16 0. 
Similarly in Fig. ^ the QMCI and II models underestimate the density in the interior region. 
The corresponding form factor is indicated in Fig. || For the diffraction radius, which 
indicates the location of the first zero of the form factor, the predictions are in reasonable 
agreement. At larger momenta a substantial model dependence in the diffraction pattern 
arises. 
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Binding energies and rms charge radii for 16 0, 4 ° Ca, 48 Ca, 90 Zr and 208 Pb are shown in 
Table [D] for various values of rj and gf. Also included are the QHD, QMCI, and QMCII 
results and the experimental values. To make a more realistic comparison with the experi- 
mental data the last row of Table [Ij indicates cm. corrections taken from Ref. []33 |. Overall 



the models QMCII, MQMC and QHD give a realistic description of the binding energies 
and radii. Including the cm. corrections the QHD and MQMC model reproduce the exper- 
imental binding energies within an accuracy of 3%. For MQMC we observe a small model 
dependence. Due to changes in the surface systematic the binding energies increase with 
rj and gl [O. The rms charge radii are insensitive to the parametrization. The QMCII 



model systematically overestimates the binding energies of the light nuclei. The error is 
maximal for Oxygen (~ 6%) and decreases with increasing mass number. We attribute this 
overbinding to the small surface energy. The binding energies in the QMCI model are too 
small. The error of the theoretical predictions ranges between 6% and 9%. 

For QMCI we could not reproduce the results for the binding energies given in Ref. || 
(see their Table 4), particularly for the light nuclei. We believe that the discrepancy is due 
to their approximate treatment of the nucleon mass. The authors of Ref. || employed a 
simple parametrization of the form 



dM* 

M*(0) « M + —— 



+ c0 2 . (35) 

4>=o 

Although the parametrization is quite accurate we found that the binding energies are 
sensitive to small variations of the fit parameter c. 

The impact of the gradient coupling is less drastic for the binding energies than for the 



nuclear shapes. To obtain a quantitative estimate of the surface energy we followed Ref. |3T] 



and fitted the error of the calculated energies 5e N (including the cm. correction) to a form 

P 



5e N = a + 



A 1 / 3 ' 



The coefficient (3 which indicates the correction to the surface energy is negative for all the 
models. We found (3 = -2.7 MeV for the QMCII model, (3 = (-1) - (-1.5) MeV for the 
MQMC models and (3 = -0.8 for the QHD model. Although the effect is small, the QMCII 
model clearly produces the smallest surface energies. 

As stated in earlier references [pHTTH the QMCI model leads to a fair description of the 
bulk properties of nuclei but gives only a poor reproduction of single-particle spectra and 
spin-orbit splittings. Spin-orbit splittings are highly correlated to the effective nucleon mass 
which is too high in the model QMCI. In view of these shortcomings the analysis of single- 
particle spectra and spin-orbit splittings is very important. The single-particle levels for 40 Ca 



are shown in Fig. [10]. The MQMC model and QHD clearly give a more realistic description 
of the energy levels than the two versions of QMC. The energies of the deeply bound states 
are too small in QMCI. Both QMCI and QMCII predict an incorrect level ordering of the 
2si/2 and ld.3/2 states (see also Refs. |],|11|]). We observe a very good agreement between 
MQMC and QHD. Generally, there is only a very weak dependence of the energy levels on 



the parameters rj and g q s JT3] . 

Results for other nuclei are similar. Spin-orbit splittings for the highest occupied proton 
and neutron states in 208 Pb and 16 O are shown in Table [Hj and Table [IV] respectively. The 
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results demonstrate directly the importance of the effective nucleon mass. The QMCI model 
systematically underpredicts the splittings [[|. The small value of the effective mass in the 
MQMC and QMCII model significantly corrects this shortcoming. 

The spin-orbit potential for 40 Ca can be studied in Fig. |Tl|. It arises in a single particle 
hamiltonian that acts on two-component wave functions. Here we follow Reinhard |M] who 
has proposed an expansion in terms of a small nucleon velocity which converges better than 
the usual Foldy-Wouthuysen reduction. In this framework the spin-orbit part is given by 



h s .o. = V s . .(r) a ■ L 
where M is defined as 

M 




dr dr 



a 



(36) 



M-UgoQ+gyVo) . (37) 



The spin-orbit potential of QMCI is rather weak explaining the small splittings in Table [D I 
and [TV]. The results for QHD and MQMC are almost identical. For QMCII the size of the 
potential is similar at the surface of the nucleus but it has only half the strength at the 
origin. 

The main results of our analysis are summarized in Table |V|. The column headings 
denote the various models discussed in this work. The rows contain specific model features 
and representative results. The first row indicates the form of the bag constant which is a 
function of the scalar field in the model MQMC (B = B((p)) and which is held at its free 
space value for the model QMCI and QMCII (B = B ). The second row contains the ratio 
of the predicted bag radius at nuclear matter equilibrium to its free space value. The specific 
form of the scalar potential is indicated in the third row followed by the equilibrium value 
of the effective nucleon mass. 

The last three rows indicate how well the models reproduce the experimental data. 
Included is the range of error for the predicted binding energies (including cm. corrections). 
As a typical representative for the quality of the spin-orbit splittings the error of the splitting 
of the ld 5 /2 and ld 3 / 2 neutron states in 40 Ca is shown. The last row indicates the integrated 
error of the charge density in 16 O. 



VI. SUMMARY 

In this paper we study properties of nuclear matter and finite nuclei based on three ver- 
sions of the quark-meson coupling model. This model describes nucleons as nonoverlapping 
MIT bags interacting through scalar and vector mean fields. The two versions QMCI and 
II differ from the third version denoted by MQMC in the treatment of the bag constant. In 
the QMCI and QMCII model the bag constant is held fixed at its free space value whereas 
in the MQMC model we assume it depends on the density of the nuclear environment. We 
employ a model for the bag constant in which the density dependence is parametrized in 
terms of the scalar mean field. The model QMCII is an extension of QMCI which includes 
additional cubic and quartic self-interactions for the scalar meson. 

The QMCII and MQMC model give rise to a sufficient number of model parameters which 
can be fit to properties of nuclear matter near equilibrium that are known to be characteristic 
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of the observed bulk and single-particle properties of nuclei. Due to the small number of 
parameters not all the desired nuclear matter properties can be fit in the QMCI model. 
Particularly, the very high value for effective nucleon mass at equilibrium is a prediction in 
this model. 

In the framework of the QMC model the nucleon size is represented by the bag radius. 
Uncertainties in the predictions for the nucleon size arise from the incomplete knowledge of 
the medium dependence of the bag model parameters. The most important quantity here 
is the bag constant. In the MQMC model an increasing bag constant leads to significantly 
swollen nucleons in the nuclear environment. In contrast, for the QMCI and II model the 
bag constant is held fixed and the bag radius decreases slightly with increasing density. In 
light of the strikingly different predictions for the nucleon size, the question arises whether 
the different models are all consistent with established hadronic phenomenology. Our basic 
goal is to study properties of nuclear matter and finite nuclei as an important test for the 
reliability of predictions for the changes of the internal nucleon structure. 

A direct relation between nuclear phenomenology and the quark picture arises from a 
redefinition of the scalar field. By introducing a new scalar field the QMC models can be 
cast into a QHD-type hadronic mean-field model with a general nonlinear scalar potential 
and a nonlinear coupling to the gradients of the scalar field. 

To make contact with the established hadronic framework we compare the QMC models 
with a QHD model including quartic and cubic scalar self-interactions calibrated to produce 
the same equilibrium properties of nuclear matter. 

Our basic result is that the models QMCII, MQMC and QHD, i.e. the models which 
can be calibrated to reproduce the empirical properties of nuclear matter, lead to essentially 
identical results in nuclear matter and for the binding energies, rms charge radii and spin- 
orbit splittings of finite nuclei. The MQMC and QMCII model reproduce large scalar and 
vector potentials and the typical size of the nonlinear potential. For the binding energies 
and spin-orbit splittings QMCII, MQMC and QHD lead to a realistic description of the 
experimental numbers. The only exception are the lightest nuclei which are systematically 
overbound in the QMCII model due to the small surface energy. The results are in accordance 
with numerous calculations in the hadronic framework. Experience has shown that an 
accurate reproduction of the nuclear matter properties leads to realistic results when the 
calculations are extended to finite nuclei |I9,|20,23|-P6[ . 

In contrast, the QMCI model does not reproduce all the empirical nuclear matter prop- 
erties. The scalar and vector mean fields are rather small and, as a consequence, the binding 
energies are systematically underestimated. The most prominent shortcoming of the model 
is the poor reproduction of single-particle spectra and spin-orbit splittings 1PHTT] . The cru- 
cial point is that this model cannot reproduce sufficient small values for the effective mass, 
which is strongly correlated with the spin-orbit force in relativistic mean-field models. 

The main difference between the models arises from the nonlinear coupling to the gradient 
of the scalar field. Such a coupling is not included in our version of QHD. It effects the surface 
energy and the nuclear shape leading to a more diffuse surface in the QMC models. The 
size of this coupling is a pure prediction which depends on the details of the underlying bag 
model. At normal nuclear matter density the gradient coupling is 6-7 times bigger for the 
QMCII model than for the MQMC model. The impact of the coupling is partly compensated 
by adjusting the scalar mass to reproduce the experimental value of the charge radius in 
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40 Ca. Although after the adjustment all the models predict nearly identical rms charge radii 
we find that the large value of the gradient coupling in the QMCII model leads to unrealistic 
features for the nuclear shapes. 

These results lead us to two basic conclusions. First, it is clear that to keep the models on 
a tractable level quark models for nuclear matter and finite nuclei always contain a number of 
parameters which cannot be determined on the fundamental level. Realistic predictions for 
the properties of finite nuclei require an accurate calibration of these parameters by fitting 
to properties of nuclear matter near equilibrium that are known to be characteristic of the 
observed bulk and single-particle properties of nuclei. This guarantees the reproduction of 
large scalar and vector potentials and a realistic description of spin-orbit splittings. Second, 
the QMC models are invoked to describe medium modifications for the internal structure of 
the nucleon. Unless more information on the fundamental level is available, the reliability 
of predictions and differences which arise in different models must be tested by analyzing 
nuclear matter and finite nuclei. Our results indicate shortcomings of the models which 
predict a decreasing size of the nucleon in the nuclear environment. The QMCI model 
cannot reproduce spin-orbit splittings on a satisfactory level. Although the QMCII model 
corrects this shortcoming, unrealistic features of the nuclear shapes arise. In principle, the 
quality of the predicted nuclear shapes can be improved by including additional gradient 
couplings. This, however, leads to new parameters which are more difficult to calibrate 
j25| . As a consequence, the formalism looses much of its predictive power and its simplicity. 



The MQMC model on the other hand, which leads to a significant increase of the nucleon 
size, provides a realistic description of nuclear matter and finite nuclei. At normal nuclear 
densities the concept of a density dependent bag constant is very useful. However, problems 
arise in an extrapolation to higher densities. Due to the sizable increase of the bag radius 
the individual bags start overlapping and the simple bag model is no longer applicable. 
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TABLES 



TABLE I. Equilibrium Properties of Nuclear Matter. The input row denotes the nuclear mat- 
ter properties which are used to calibrate the models MQMC, QMCII and QHD. The second row 
indicates the nuclear matter properties of the QMC model which contains only three free parame- 
ters. 







P° 


M*/M 


e 


K 




a 4 


Input 


1.3 far 1 


0.1484 far 3 


0.63 


-15.75 MeV 


224.2 MeV 


35 


MeV 


QMCI 


1.3 far 1 


0.1484 far 3 


0.80 


-15.75 MeV 


280.0 MeV 


35 


MeV 



TABLE II. Binding energy per nucleon e N (in MeV) and rms charge radius r c in (fm) for several 
closed shell nuclei. 



Model 


m s 


16Q 




40 


Ca 


48 


Ca 


90 


Zr 


208 


Pb 


V 9l 


(MeV) 


e N 


r c 


e N 


r c 




r c 


e N 


r c 


e N 


r c 


3 2 


543 


-7.26 


2.74 


-8.15 


3.48 


-8.22 


3.50 


-8.35 


4.29 


-7.57 


5.56 


3 3 


544 


-7.28 


2.74 


-8.20 


3.48 


-8.27 


3.50 


-8.38 


4.29 


-7.59 


5.56 


4 2 


546 


-7.37 


2.74 


-8.26 


3.48 


-8.33 


3.50 


-8.42 


4.29 


-7.63 


5.56 


4 3 


547 


-7.41 


2.74 


-8.30 


3.48 


-8.37 


3.50 


-8.45 


4.29 


-7.65 


5.56 


QMCI 


464.5 


-6.65 


2.75 


-7.73 


3.48 


-7.59 


3.53 


-7.85 


4.31 


-7.13 


5.57 


QMCII 


685 


-7.83 


2.74 


-8.50 


3.48 


-8.58 


3.52 


-8.61 


4.30 


-7.74 


5.56 


QHD 


500.8 


-7.18 


2.74 


-8.14 


3.48 


-8.19 


3.50 


-8.30 


4.29 


-7.57 


5.56 


Exp. 




-7.98 


2.73 


-8.45 


3.48 


-8.57 


3.47 


-8.66 


4.27 


-7.86 


5.50 


CM. 




-0.67 




-0.21 




-0.18 




-0.08 




-0.03 





TABLE III. Spin-orbit splittings of the highest occupied proton and neutron levels in Pb. 
For the MQMC model the scalar coupling to the quarks is g% = 3. 



protons 




rj = 3 


7] = 4 


QHD 


QMCI 


QMCII 


expt. 


35] 


AE(2d 5/2 - 2d 3/2 ) 


(MeV) 


-1.42 


-1.41 


-1.39 


-0.55 


-1.55 


-1.3 


AE(lg 9/2 - lg 7/2 ) 


(MeV) 


-3.40 


-3.40 


-3.43 


-1.22 


-3.21 


-4.0 


neutrons 




rj = 3 


r\ = 4 


QHD 


QMCI 


QMCII 


expt. 


m 


AE(3p 3/2 - 3pi/ 2 ) 


(MeV) 


-0.68 


-0.68 


-0.66 


-0.26 


-0.74 


-0.9 


AE(2f 7/2 - 2/ 5/2 ) 


(MeV) 


-1.80 


-1.78 


-1.74 


-0.74 


-2.03 


-1.8 
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TABLE IV. Spin-orbit splittings of the highest occupied proton and neutron levels in 16 0. For 
the MQMC model the scalar coupling to the quarks is <?f = 3. 



protons 




77 = 3 


77 = 4 


QHD 


QMCI 


QMCII 


expt. 


35] 


A£(lp 3 /2 - IP1/2) 


(MeV) 


-5.16 


-5.18 


-5.27 


-1.87 


-4.71 


-6.3 


neutrons 




r\ = 3 


■x] = 4 


QHD 


QMCI 


QMCII 


expt. 35] 


AE(lp 3/2 - lpi/2) 


(MeV) 


-5.22 


-5.24 


-5.34 


-1.88 


-4.75 


-6.1 





TABLE V. Model Summary. For a detailed description see text. 



Model 




MQMC 




QMCI 


QMCII 


QHD 


V\9l 


3\2 


3\3 


4\2 


4\3 








Bag constant 


B{<P) 


B{<P) 


B(<P) 


B{4>) 


Bo 


B 




R/Ro 


1.37 


1.28 


1.37 


1.28 


0.99 


0.97 




Seal, potential 


\m 2 s 4> 2 




\mlct) 2 




\m 2 s (t> 2 


(\ml4> 2 + jt^ 




M*/M 


0.63 


0.63 


0.63 


0.63 


0.80 


0.63 


0.63 


Binding energies 


1-3% 


1-3% 


1-3% 


1-2% 


6-9% 


1-7% 


1-3% 


Spin-orbit splitting 


12% 


12% 


12% 


12% 


68% 


18% 


11% 


Charge density 


2.6% 


2.6% 


2.5% 


2.4% 


3.0% 


3.7% 


2.5% 
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FIGURES 




FIG. 1. Bag radius as a function of the density for the models MQMC, QMCI and QMCII. 
The solution for the QMCII model terminates at /9 max ss 1.23/9°. 
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FIG. 2. Effective nucleon mass as a function of the density. 
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FIG. 3. Predicted nonlinear scalar potential as a function of the transformed scalar field 
g $ = M-M*. 
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FIG. 5. Gradient coupling as a function of the transformed scalar field g & = M — M*. The 
function is identical for the models QMCI and QMCII. 
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(a) (b) 



FIG. 6. Charge density for 40 Ca. In part (a) the same mass for the scalar meson 
(m s = 500.8MeV) was used for all the models. In part (b) the scalar mass was adjusted to 
reproduce the experimental charge radius in 40 Ca. The parameters for the MQMC model are 
rj = 4 and g q s = 3. 
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FIG. 7. Isotope shift in calcium. The difference between the Ca and Ca charge densities is 
shown. The experimental (shaded) curve is taken from Ref. |p^| . The parameters for the MQMC 
model are r\ = 4 and gf = 3. 
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FIG. 8. Charge density for 16 0. The parameters for the MQMC model are r\ = 4 and gj = 3. 
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FIG. 9. Form factor for 16 0. The parameters for the MQMC model are r] = 4 and g q s = 3. 
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FIG. 10. Single particle spectrum of 40 Ca. The parameters for the MQMC model are rj 
and gl = 3. The experimental numbers are taken from Ref. 
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FIG. 11. Spin-orbit potential for 40 Ca. The parameters for the MQMC model are rj = 4 and 
= 3. 
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